Nuclear magnetic resonance as a probe of electronic states of Bi 2 Se 3 
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We present magnetotransport and 209 Bi nuclear magnetic resonance (NMR) data on a series of 
single crystals of Bi2Se3, Bi2Te2Se and Cu a; Bi2Se3 with varying carrier concentrations. The Knight 
shift of the bulk nuclei is strongly correlated with the carrier concentration via a hyperfine coupling 
of 27 /ieV, which may have important consequences for scattering of the protected surface states. 
Surprisingly we find that the NMR linewidths and the spin lattice relaxation rate appear to be 
dominated by the presence of localized spins, which may be related to the presence of Se vacancies. 

PACS numbers: 76.60.-k, 72.80.-r, 31.30.Gs 



I. INTRODUCTION 

Topological insulators are a novel form of condensed 
matter in which the bulk is electrically insulating but the 
surface electronic states remain gapless and conducting. 1 
This unusual situation can emerge in materials with 
large spin-orbit couplings such that the topology of the 
bulk band structure differs from that of the surrounding 
material. 2 At the boundaries of the topological material 
the surface electron states exhibit a Dirac dispersion, and 
the number of such Dirac points in k-space is odd. Fur- 
thermore, the strong spin orbit coupling locks the elec- 
tron spin direction to the momentum giving rise to chiral 
states. 3 These chiral surface states are protected, mean- 
ing that perturbations such as defects and disorder on the 
surface will not alter their character. More importantly 
electron backscattering is strongly suppressed, and may 
enable dissipationlcss polarized spin currents which could 
prove useful in spintronics applications. 4 

In 2009 the first of the bismuth chalcogenide family of 
materials, Bi2Se3, was discovered to exhibit surface states 
with the single Dirac cone characteristic of a 3D topolog- 
ical insulator. 5 These compounds have attracted signifi- 
cant attention because large single crystals can be grown 
and clean surfaces can be exposed by cleaving, which 
enables detailed studies of the electronic dispersion via 
angle resolved photocmission spectroscopy (ARPES). 1 
Bi2Se3 has a rhombohcdral crystal structure consisting 
of hexagonal stacked planes of Bi and Se (see Fig. 1). 
The planes are organized into quintuple layers, and each 
layer has Se atoms on top and bottom. In principle this 
material should be an intrinsic semiconductor with a gap 
of 350 meV; however, there is an inherent chemical ten- 
dency for Se atoms to be missing from the lattice. 6 These 
vacancy sites act as donors, raising the chemical poten- 
tial out of the insulating regime and forming a degen- 
erate semiconductor. In nearly all cases the vacancies 
place the Fermi level well in the conduction band, result- 
ing in metallic behavior. As a result the exotic trans- 
port properties of the surface states are often masked by 
the parallel bulk conduction channel. 7 Although several 
techniques have emerged to dope these materials so that 
the bulk states become insulating and the surface states 



are tuned to the Dirac point, relatively little is known 
about the microscopic electronic response of both the 
bulk and surface states. 8 ~ 10 In order to investigate the 
bulk states in more detail we have conducted 209 Bi nu- 
clear magnetic resonance (NMR) and magnetotransport 
studies of Bi 2 Sc 3 grown in excess Se, Bi 2 Tc 2 Sc, and Cu 
doped Cu.j;Bi2Sc3. We find clear evidence of correlations 
between the NMR response and the carrier concentration 
in the Bi2Sc3 samples. 

Despite extensive research on the bismuth chalco- 
genides, there have been relatively few studies of the 
NMR response of this material. 209 Bi spectra and 
relaxation rates were reported in bulk Bi2Sc3, n and 
125 Te NMR has been reported in nanoscale powders of 
Bi2Tc3. 12 Depending on the surface to volume ratio, nu- 
clei at the surface may have sufficient spectral weight 
to contribute to the NMR signal, and in the former 
study the spectra revealed features consistent with sur- 
face nuclei. 13 Moreover, the spin-lattice relaxation rate 
(T-f x ) of these surface nuclei was enhanced compared to 
the bulk. The authors concluded that this enhancement 
was consistent with metallic protected surface states. 
The spin-lattice relaxation rate of the nuclei may be 
driven by spin-flip scattering between the electron and 
nuclear spins via the hyperfine interaction. However, it 
is unclear whether T-f" is enhanced or suppressed by the 
exotic surface states, or whether the hyperfine interaction 
would have sufficient magnitude. In order to understand 
the interaction of nuclei with the topological states it is 
essential to understand the relaxation rates and hyperfine 
couplings of both the bulk and surface states. 

The bismuth chalcogenides are well suited for NMR be- 
cause Bi, Se and Te all have NMR active isotopes: 209 Bi 
(I = §, 100% abundant), 77 Se (7 = ±, 7.63% abundant) 
and 125 Tc (I = |, 7.07% abundant). 14 209 Bi is particu- 
larly useful because not only is the signal intensity larger 
than that of 77 Se, but it is also sensitive to the local 
charge environment via the quadrupolar interaction. The 
nuclear spin Hamiltonian for 209 Bi is given by: 
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where 7 = 0.6842 kHz/G is the gyromagnetic ratio, 
Hq = 9 T is the external field (except for sample #1 and 
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FIG. 1. A representative powder X-ray diffraction pattern 
of Bi2Se3 (black line) compared with theoretical prediction 
(blue) Inset: Crystal structure of Bi2Se3. Bi (blue) and Se 
(green and yellow) atoms form quintuple layers. 



the Bi 2 Te 2 Se in which case H = 11.724 T), I a are the 
nuclear spin operators, v cc is the component of the elec- 
tric field gradient (EFG) tensor along the c-direction, 77 
is the asymmetry parameter of the EFG tensor, and T-Lm 
is the hypcrfinc interaction between the Bi nuclear spins 
and the electron spins. 15 Since there is a single Bi site in 
the unit cell with axial symmetry, r\ is zero, and the spec- 
trum consists of a set of nine equally spaced resonances 
separated by v cc . The hyperfine coupling can be written 
as %hf = 7%/iBj4I • S, where S is the electron spin (of 
either the bulk carriers or of the surface state electrons) 
with g-factor g, and A is the direct contact hyperfine in- 
teraction. The hyperfine coupling gives rise to Knight 
shift K = K + Ax, where \ is the bulk magnetic sus- 
ceptibility and Kq is the temperature independent orbital 
shift. For a degenerate semiconductor or metallic system, 
the susceptibility is given by \ = ff 2 Ms-^(0), where N(0) 
is the density of states at the Fermi level. The nuclear 
spin Hamiltonian (Eq. 1) gives rise to nine resonances of 
the Bi nuclei in Bi2Se3 at frequencies: 



1 H {1 + K) 



(2) 



where n 
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II. EXPERIMENTAL METHODS 

Several different single crystals of Bi2Sc3 were prepared 
by the Bridgman method from varying mixtures of el- 
emental Bi and Se in evacuated quartz ampoules. 16 A 
vertical Bridgman furnace was set up to lower an am- 
poule through a temperature gradient of about 9 °C/cm 
around the melting point of 710 °C at a rate of 2 mm/hr. 
Stoichiometric mixtures of the elements were pre-mcltcd 
and homogenized in the ampoules at 800 °C for 12 hours 
and then furnace cooled. The ampoule was then placed 



into the Bridgman furnace and lowered through the gra- 
dient. Samples that were prepared with the initial mix- 
ture having a nominal stoichiometry of Bi2Se3 formed an 
ingot with large columnar single crystals of a preferred 
orientation. However, samples grown with an excess of 
Se formed polycrystalline ingots, from which only a few 
single crystal samples could be obtained. We found that 
samples grown simply by annealing in quartz ampoules 
facilitated the formation of large single crystals, though 
not with a preferred orientation. The Cu-dopcd sam- 
ple grown following the method described in Reference 
9. Powder X-ray diffraction measurements confirmed the 
phase purity of each sample (see Fig. 1). Table I sum- 
marizes the growth methods and properties of all of the 
samples investigated. 

A low speed diamond wheel saw was used to cut rect- 
angular bar shapes from the large single crystals in our 
ingots of Bi2Se3. Standard 4- wire resistivity measure- 
ments and 5-wire Hall measurements were performed us- 
ing a Quantum Design Physical Properties Measurement 
System (PPMS), and also with a Keithley current source 
and nanovoltmetcr in a closed-cycle refrigerator. Magnc- 
toresistance measurements were performed in the PPMS 
by using the 4-wire method with silver painted contacts 
to measure sample resistance at a temperature of 1.9 K. 
Shubnikov-dc-Haas oscillations were observed in some, 
but not all, samples. The carrier concentrations and mo- 
bilities were estimated from the Hall coefficients and re- 
sistivities. The NMR spectra and relaxation rates were 
measured using a standard spin echo pulse sequence at 
varying frequencies in fields of 9 T and 11.7 T. 



III. RESULTS AND DISCUSSION 
A. Charge Transport 

The electronic properties of our samples were strongly 
dependent on growth conditions. The Hall voltage and 
resistivity of several samples were measured as a function 
of temperature, and these results are shown in Figs. 2a 
and 2b. The carrier concentration n was estimated from 
the Hall coefficient, Rh using the formula Rh = —1/nq, 
where q = — e is the charge of the electron carriers (from 
the Se vacancies) . The Hall coefficients of samples from 
batches B and E have carrier concentrations of order 10 19 
cm" 3 . Samples grown with excess Se (batches C and D) 
have lower carrier concentrations on the order of 10 18 
cm -3 . Excess Se can reduce the carrier concentration by 
reducing the number of Se vacancies, which donate two 
electrons each. On the other hand crystals grown in ex- 
cess Se can also have Bi/Se antisite defects, which donate 
a single electron. 17 Although our NMR spectra indicate 
that we introduce some disorder in the samples with ex- 
cess Se, it is not readily apparent that this disorder arises 
from antisite defects. 

Figures 3 and 4 display the carrier concentration and 
lattice constants versus the nominal percentage of Se in 



TABLE I. Measured Bi2Se3 samples, and their growth methods and properties at 10 K. The atomic percentage of Se for a 
stoichiometric Bi2Se3 mixture is 60%. 



Sample Batch Growth mechanism % Se n (cm ) /i(V"m 



/sdH (T) 



K (%) Ac^ ctr (kHz) Ai^Q (kHz) 



(kHz) 



#1 


A 


Ampoule anneal 


60 


N/A 


N/A 


N/A 


0.647 ± 0.002 


38 ±3 


25 ±3 


158.3 ±0.6 


#2 


B 


Bridgman 


60 


1.58 x 10 19 


0.540 


127 ±6 


0.619 ±0.001 


30 ±2 


30 ±2 


160.5 ±0.5 


#3 


C 


Bridgman 


65 


N/A 


N/A 


76 ±4 


0.51 ±0.01 


86 ±4 


N/A 


168 ±2 


#4 


c 


Bridgman 


65 


6.38 x 10 17 


1.049 


38 ±4 


0.356 ± 0.003 


80 ± 10 


N/A 


159.9 ±0.06 


#5 


D 


Bridgman 


62.5 


3.46 x 10 18 


0.249 


N/A 


0.441 ± 0.004 


56 ±5 


11 ± 16 


166 ±2 


#6 


E 


Bridgman 


6(1 


1.35 x 10 19 


0.108 


N/A 


0.667 ±0.001 


37 ±2 


43 ± 1 


150.5 ±0.5 


#7 


B 


Bridgman 


6(1 


9.41 x 10 18 


0.645 


N/A 


N/A 


N/A 






#8 


B 


Bridgman 


6(1 


N/A 


N/A 


N/A 


N/A 


N/A 






#9 


B 


Bridgman 


60 


1.91 x 10 19 


0.680 


N/A 


N/A 


N/A 
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FIG. 2. (a) Carrier concentrations of Bi2Se3 as determined 
by Hall measurements. The stoichiometric samples (#2 and 
#6 marked as in Fig. 7) show the highest concentrations, 
as well as #7 and #9 (hollow diamonds and slanted bars, 
respectively) from batch B. Excess Se (samples marked as in 
Fig. 7) can lower the concentration by an order of magnitude, 
(b) Temperature dependence of resistivities of our samples. 
Sample =#=4 results are from the closed-cycle refrigerator; all 
other samples are from the measurements in the PPMS. 



these samples. It is clear that adding Se reduces the num- 
ber of electrons, consistent with expectation. Further- 
more, the monotonic variation of the lattice constants is 
consistent with Vegard's Law, and indicates that the Se 
is doping homogeneously in the materials. Curiously, the 
behavior of Bi2Se3 contrasts with that of Bi2Te3 (Fig. 3), 
in which the carrier concentration first decreases with in- 
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FIG. 3. Carrier concentration of Bi2Se3 (•) and Bi2Te3 (o) 
as a function of the atomic percentage of Se/Te in the initial 
mixture of Bi and Se. Data for the Bi2Te3 reproduced from 
18. 



creasing Tc concentration and then increases. 18 The dif- 
ference is due to the different types of defects that occur 
in the two materials: in stoichiometric Bi2Se3 the domi- 
nant defect is Se vacancies, so the addition of Se reduces 
the vacancies and may even add Bi vacancies which act as 
acceptors. In stoichiometric Bi2Te3 the dominant defect 
is a substitution of Bi for Te which acts as an acceptor; 
excess Te causes the dominant defect to be substitution 
of Te for Bi, a donor. 

The resistivity, p, of these samples revealed metallic be- 
havior, consistent with a high, temperature independent 
carrier concentration (Fig. 2b) as expected for a degener- 
ate semiconductor. We extract the mobilities, y, = 1/pne 
from the resistivity and Hall measurements, as shown in 
Fig. 5. These values are between 1-10 m 2 /V-s and are 
comparable to those in the literature. 7 



B. Quantum Oscillations 

The properties of the doped carriers can be further 
investigated by magnetoresistance measurements. The 
resistivity at 1.9 K was measured up to a field H = 9 




28.63 



60 61 62 63 64 65 
Nominal Atomic Percent Se 



FIG. 4. Lattice constants a and c as a function of the atomic 
percentage of Se. Solid lines are best linear fits to the data. 
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FIG. 5. The mobilities versus temperature as determined 
from Hall constant and resistivity measurements. 



T along the c-axis and perpendicular to the applied cur- 
rent, and the data were fit to a polynomial of degree 
five. The difference, Ap, between the measured values 
and the polynomial fit display clear Shubnikov-dc Haas 
(SdH) oscillations, as shown in Fig. 6. The oscillations 
are described by the Lifshitz-Kosevitch formula: 



Ap = po sin 



F 

27r ( Jj - 75 



(3) 



where H is the magnetic field, F is the SdH frequency 
related to the Fermi surface area, and jg is a phase factor 
related to the Berry phase. 19 The amplitude po depends 
on the temperature, the scattering rate, and the Dingle 
temperature (defined as ft/27rfcsT, where r is the electron 
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FIG. 6. Ap versus inverse field, H^ 1 at 1.9 K in Bi2Se3, 
where Ap is the difference between the measured resistivity 
and a polynomial fit as described in the text. The curves are 
offset vertically for clarity. NMR was performed on the sam- 
ples which showed oscillations (marked as in Fig. 7). Other 
samples from batch B (samples #7, #8 and #9) showed no 
oscillations. 



scattering time). From the data we estimate the Fermi 
wavevectors using the formula: 



2?re 



(4) 



where A is the maximal cross-sectional area, and assum- 
ing the Fermi surface to be spherical. Only one sample 
from batch B showed such oscillations, with a frequency 
of about 124 T. Two other samples grown with excess Se 
revealed lower SdH frequencies of 72 T and 35 T (Fig. 
6). These lower frequencies are consistent with the Fermi 
surface having a smaller maximal cross-section due to a 
lower concentration of electron donors. These oscillations 
correspond to Fermi wavevectors kp ~ 0.02 — 0.06A -1 . 
These values are roughly the same order of magnitude 
as those reported by ARPES measurements of samples 
grown by Se self-flux (0.04 A -1 ) and traveling floating 
zone solvent growth (0.08 A" 1 ). 11 Although we do not 
have data on the temperature dependence of the SdH 
amplitudes, we can estimate the upper limit of the effec- 
tive mass at 1-4 times the bare electron mass based on 
a fit of the pcak-to-pcak amplitudes of each cycle versus 
the inverse field at the midpoints of the cycles, and the 
measured values of the mobility. 



C. Nuclear Magnetic Resonance 

The Bi NMR spectra were measured at 10 K in several 
different single crystal samples and are summarized in 



Fig. 7. The spectra were fit to a sum of nine Lorentzians: 

4 



i. — / 7]- 



Aw n 



(w - u>„) 2 + Aw 2 



(5) 



where the frequencies u;„ are given by Eq. 2, the 
linewidths are given by Aw n = -^(^a/) 2 + n((5wg) 2 , 
and Slum and <5u;q are the magnetic and quadrupolar con- 
tributions. This equation accounts for both a magnetic 
and quadrupolar broadening for each transition. 21 The 
amplitudes A n were allowed to float, although in princi- 
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FIG. 7. 209 Bi NMR spectra of Bi 2 Se 3 single crystals in a field 
of 9 T with Ho || c at 10 K (except for #2 at 20 K). The 
quadrupolar splitting appears to be sample-independent, but 
the Knight shift depends on carrier concentration. Samples 
#1 (•); #2 (■), and #6 (EB) all show a similar spectral profile 
and were all made from stoichiometric mixtures. Samples 
#3 (a) and #4 (+) have pronounced outer peaks, consistent 
with observations in samples grown in Se self-flux. 11 Sample 
#5 (x) shows a more typical profile. Disorder broadens the 
quadrupolar satellites and washes out the spectra in the 281 
K Bi 2 Te 2 Se (V) and Cu^BigSea (♦) samples. 20 



pie they are coupled and related to the spin echo decay 
rate. 11 The fitted values of K, Sojm, Sujq and v cc are 
reported in Table 1. For all but Bi 2 Te2Se, the nine satel- 
lites of the 209 Bi are distinguishable. The EFG splitting 
v C c ^150 — 170 kHz is somewhat smaller than previous 
experiments on other bismuth compounds. 22 ' 23 For some 
of the spectra, the intensity of the satellites follow an 
atypical distribution in which the central resonances are 
suppressed relative to the satellites. This effect has been 
observed previously in Bi 2 Se3 and explained in terms of 
a spin-spin decoherence rate that depends on the partic- 
ular nuclear transition. 11 Not all of the spectra exhibit 
this unusual distribution, however. In particular, one of 
the samples grown with excess Se shows a more typical 
spectral profile. This effect probably reflects a longer T 2 
relaxation time for these samples. Aside from sample 
#5, this feature appears to be correlated with the carrier 
concentration and is more pronounced at lower n. 

Figure 8 shows the Knight shift, K, versus the carrier 
concentration, n, as determined from the Hall constant. 
The Knight shift arises because spins of the electron car- 
riers couple to the nuclei via the hyperfine interaction. 
In order to estimate the magnitude of this coupling, we 
fit the data to the the expression K = K + A\ m , where 
the molar susceptibility is given by the Pauli expression: 



Xn 



-ir 2 



h 2 TT 2 



(37T 2 n) 



1/3 



(6) 



Here n is the carrier concentration as determined by the 
Hall constant, m* = 0.19mo is the effective mass as de- 
terminted by ARPES measurements (too is the bare elec- 
tron mass), V m = 258 cm 3 /mol is the molar volume, 
and A is the hyperfine coupling. 24 This expression fits 
the data well (solid line in Fig. 8) with fit parameters 
K = 0.19 ± 0.06%, and A = 9.7 ± 1.6 MOe/^s- In 
units of energy, 209 jhA^B = 27/ieV; or alternatively a 
hyperfine field of Bhf ~ 735 G experienced by the elec- 
trons. This value is larger than typical transferred hyper- 
fine couplings, but is comparable to other semiconductors 
with on-site Fermi contact couplings. 25 ' 26 It is likely that 
the hyperfine interaction for the surface states would be 
similar to that of the bulk. The large value of the hy- 
perfine coupling as well as the high natural abundance of 
209 Bi may limit the coherence time of any spin polarized 
currents, but it is unclear how the spin lattice relaxation 
process would be modified because of the protected na- 
ture of the surface states. 27 

Using this hyperfine constant we can estimate the 
carrier concentration of the Cu doped sample at n ~ 
2.5 x 10 19 cm -3 for our sample with a nominal Cu con- 
tent of a; = 0.12. The carrier concentration remains an 
order of magnitude lower than that of superconducting 
samples, and probably explains the absence of supercon- 
ductivity in these samples. 9 Based on the spectrum seen 
in Fig. 7 it is clear that there is broad distribution of 
local EFGs that affect the quadrupolar satellites. This 
disorder may result from either intercalation of Cu be- 
tween Bi 2 Sc3 layers or via Cu/Bi substitution. 
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FIG. 8. Knight shift, K, versus carrier concentration, n. The 
positive correlation is consistent with a hyperfine coupling to 
electrons. The solid line is a fit as described in the text. 
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FIG. 9. The magnetic contribution to the linewidth, SoJm 
versus carrier concentration (lower axis) and versus Se con- 
centration (upper) axis. 



Interestingly, there appears to be a negative correlation 
between the carrier concentration and the magnetic con- 
tribution to the linewidth, 5ujm, as shown in Figure 9. 
The intrinsic linewidth of the Bi should be determined 
by the second moment of the nuclear dipole-dipole inter- 
action, which is on the order of 500 Hz, two orders of 
magnitude lower than our observations. 15 Defects such 
as impurities or vacancies in the crystal can increase this 
linewidth by creating a local variation in the chemical 
shift, Kq, as well as giving rise to local magnetic fields 
from localized moments. The fact that that the atypical 
lineshapes appears to correlate with carrier concentra- 
tion further suggests the presence of magnetic impurities, 
since fluctuating magnetic fields enhance the spin deco- 
herence rate, T 2 ~ , and suppress the signal intensity of 
the central lines. Although excess Se would be expected 
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FIG. 10. Spin-lattice relaxation rates of the 209 Bi versus tem- 
perature for various samples with different carrier concentra- 
tions. The open data points are reproduced from Ref. 11, and 
exhibit a strong temperature dependence. The solid lines are 
predictions based on the Korringa formula as described in the 
text. 



to narrow the lines by reducing the Se vacancies, in our 
experiment adding Se seems to broaden the spectra. This 
result suggests that excess Se may be introducing lattice 
defects other than Se site vacancies, which may be giv- 
ing rise to localized unpaired electrons. Excess Se may 
occupy positions normally occupied by Bi as discussed 
previously, or it may cause Bi vacancies to form which 
act as acceptors, or it may occupy interstitial positions 
in the lattice. 28 On the other hand, the broadening of the 
linewidth for low n may reflect an inhomogeneous distri- 
bution of local carrier concentrations on a macroscopic 
scale. This inhomogeneity may be more pronounced at 
small n; however this scenario does not explain the atyp- 
ical lineshapes. 

The spin-lattice relaxation rate was measured at the 
central transition of several of the samples, and the mag- 
netization recovery was fit to the standard expression for 
spin 9/2 nuclei: 

Mit) = Mo [1 - 2/ (-^Le^stm + i56| e - 2 8*/T 1 + 

J5_p-15t/Ti , 24 -6t/Ti , J_ p t/T 1 \] (n\ 

65 e ' 715 e ^ 165 e )\ ' V l 

where the equilibrium magnetization, Ma, the inversion 
fraction, /, and T\ are fitting parameters. 29 The results 
are shown in Fig. 10. Although our observations are sim- 
ilar in magnitude to other published data, we find a much 
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FIG. 11. Recovery of the Bi echo after inversion recovery at 
the central line in sample #2 at 10 K. The solid line is a fit 
to Eq. 7 and the dotted line is a fit to Eq. 9. 



weaker temperature dependence. 11 In principle there are 
three possible mechanisms for spin lattice relaxation in 
a nondegenerate semiconductor: (a) Korringa relaxation 
via the contact hyperfine interaction with the conduc- 
tion electrons, (b) spin diffusion from localized electron 
spins, and (c) quadrupolar relaxation via phonons. We 
first consider case (a), in which case the spin- lattice re- 
laxation rate should be given by: 



t; 



(K - K ) 2 T/ K 



(8) 



where the Korringa constant is given by k = 
.g 2 ^/47rfcs7 2 and g± = 23 is the g-factor for the elec- 
tron carriers in this material for the field perpendicular 
to the c-axis. 30 Note that although the field lies in the c- 
direction, the fluctuating hyperfine field driving the spin- 
flip scattering lies in the plane. The calculated values of 
T[~ arc shown as solid lines in Fig. 10, and are about 
an order of magnitude lower than the measured values. 
It is possible that the discrepancy can be attributed to 
an overestimate of the orbital shift Kq from the fit to 
Eq. 6 in Fig. 8. In fact, samples #3, #4 and #5 exhibit 
linear behavior above 60 K, as expected for Korringa re- 
laxation. On the other hand, for samples #1, #2, and 
#6 the relaxation the temperature dependence is sublin- 
ear, suggesting that the Korringa mechanism is not the 
dominant relaxation channel for the samples. 

Case (b), spin diffusion from impurities, is a possi- 
bility because 209 Bi is 100% abundant, and Sc vacan- 
cies/interstitials can bind localized electron spins that are 
potential sources of scattering/relaxation. In this case 
the magnetization recovery typically exhibits a distribu- 
tion of relaxation rates because spin diffusion drives the 
relaxation of Bi nuclei with varying distances. Figure 11 
shows the recovery of the magnetization for sample #2 



M = M 



1 - 2/e 



<t/TxY 



(9) 



where the stretched exponent P is a measure of the width 
of the Tj -1 distribution. 31 The data fit better to the 
stretched exponential form. The temperature depen- 
dence of T^ 1 using Eq. 9 is essentially identical to that 
of the conventional magnetic relaxation, Eq. 7, except 
that Tj" 1 values are about a factor of 50 times larger. 
However, it is important to note that another explana- 
tion for the poor fit to Eq. 7 is that multiple transitions 
may have been inverted by the pulse sequence, rather 
than solely the central transition. As seen in Fig. 7, the 
quadrupolar splitting is small (~ 150 kHz), and the exci- 
tation bandwidth of 167 kHz could have partially excited 
some of the satellite transitions, especially if one consid- 
ers that the line width is comparable to the splitting. 
As a result the initial conditions of the magnetization 
recovery would be affected and hence the recovery func- 
tion (Eq. 7) would be modified. Furthermore, mutual 
spin flips between Bi neighbors, which drives spin diffu- 
sion, are likely suppressed because of the finite quadrupo- 
lar splitting. 32 The bulk magnetic susceptibility of these 
samples is diamagnetic, with a clear Curie tail for tem- 
peratures below 50 K. The local moments responsible for 
this Curie contribution may also be the source of the en- 
hanced magnetic linewidth seen in Fig. 9 as well as the 
stretched behavior of the relaxation. However there is no 
evident correlation between T\ , Sum and the Curie con- 
stant, suggesting that either there are extrinsic phases 
contributing to the susceptibility or the relationship be- 
tween these quantities is complex. Therefore it is unclear 
which of these two scenarios is the correct explanation for 
the stretched nature of the magnetization recovery. 

Case (c), quadrupolar relaxation, is possible because of 
the large quadrupolar moment of the Bi. This mechanism 
is often difficult to discern, and may occur in combination 
with magnetic mechanisms. 33 In this case, the tempera- 
ture dependence of T[~ is driven by changes in the lattice 
and depends on the phonon spectrum. One might expect, 
then, that T^ would correlate with the presence of any 
lattice imperfections. However, T^ docs not appear to 
depend on the electronic mobility, suggesting that this 
mechanism is not dominant. 



IV. CONCLUSIONS 

NMR and transport measurements of a series of Bi2Sc3 
crystals reveal a strong variation in the Knight shift and 
linewidth as a function of carrier concentration. By an- 
alyzing the Knight shift we determine a hyperfine cou- 
pling constant A ~ 27/j.eV. This coupling leads to a large 
Knight shift in these samples with high carrier concentra- 
tions. Because this contact hyperfine interaction is de- 
termined by the magnitude of the electron wavefunction 



at the nucleus, it is strongly controlled by large energy 
scale atomic parameters. Therefore it should not vary 
significantly from one chalcogenide to another. Further- 
more it should not be significantly different for Bi located 
on the surface rather than in the bulk. As a result this 
electron-nuclear interaction can give rise to scattering of 
the protected surface state electrons and possibly play 
a role in dephasing spin polarized currents. 34 The spin- 
lattice relaxation rate we observe is surprisingly large 
and weakly temperature dependent. Furthermore, the 
magnetic broadening we observe increases with decreas- 
ing carrier concentration. Our results may suggest the 
presence of localized spins, possibly associated with Se 



vacancies or interstitial donor sites that are partially ion- 
ized. 
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